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We study theoretically competition between the charge ordering and the superconductiv- 
ity in the two-dimensional organic conductor /3-(77ieso-DMBEDT-TTF)2PFg. We analyze the 
extended Hubbard model on a weakly-dimerized lattice based on the random phase approx- 
imation and Eliashberg equations. We found the reentrant behavior of the checkerbord-type 
charge-ordered phase in the phase diagram, and the triplet superconductivity due to the 
charge fluctuation in the neighboring region. In low temperatures, the singlet superconduc- 
tivity also appears due to the enhancement of the spin fluctuation. 
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1. Introduction 

Organic conductors have attracted attention for many years. About two decades ago, 
superconductivity was found in the low dimensional organic conductor (TMTSF)2PF6 at the 
critical temperature ~ 0.9 K under the pressure P ~ 12 kbar. Since then, about 130 
kinds of organic superconductors have been discovered. It is notable that BEDT-TTF salts 
(ET salts) are about 50 kinds of them.^'^ The ET salts are mainly quasi-two-dimensional 
(2D) conductors with 3/4-filled vr band of the donor molecules. These materials exhibit very 
interesting electronic properties such as superconductivity, magnetism and charge ordering. 

For example, k-(ET)2X salts show superconducting (SC) transitions next to the anti- 
ferromagnetic state in the (P, T) phase diagram. The theoretical calculations based on the 
half-filled extended Hubbard model where the strong dimerization is assumed show that this 
superconductivity is mediated by antiferromagnetic spin fluctuation.'^"^ In a-(ET)2l3 salt, 
there are several theoretical works for the superconductivity in the presence of a charge order- 
ing by the random phase approximation (RPA).'^''^ In this case, the SC state is also considered 
to be mediated by the antiferromagnetic spin fluctuation. For 0-(ET)2X salts where the var- 
ious charge-ordered (CO) states are observed, there are several theoretical works on the CO 
states by the mean field theory.^' ^'^ The superconductivity of this material is explained that 
it is induced by both the spin and the charge fluctuations.^^ 
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In this paper, wc focus on /3-(TOes(>DMBEDT-TTF)2PF6 salt (/?-(DMeET)2PF6 salt) 
which has been recently synthesized and investigated by Kimura and Mori et al}'^^^'^ The 
conduction layer of this material consists of weakly dimerized two molecules per a unit cell 
with the 3/4- filled band. This salt shows the metal- insulator (MI) transition at Tmi ~ 90 K 
without the anomaly of the magnetic susceptibility around Tmi- The insulating state exhibits 
the checkerboard- type charge ordering (see Fig. 1(b)) where the charge disproportion occurs 
within the dimer unit. Estimating the nearest neighbor Coulomb interaction by the point 
charge approximation, it is reported that the checkerboard type charge ordering cannot be 
explained, so that the effective molecular size and the electron-phonon interaction should be 
taken into account. The most remarkable point is that the SC state of /3-(DMeET)2PF6 is 
next to the CO state in the (P, T) phase diagram, with the critical temperature Tc ~ 4 K 
under the pressure P ~ 4 kbar.^^'^^ Tc decreases with increasing the magnetic field, and the 
superconducting transition is completely suppressed above 4 T down to 0.5 K. 

It is well known that SC transitions often appear next to the magnetic ordered region of 
the {P,T) phase diagram not only in organics but also in heavy fermion systems. This 
has been considered as a typical character of superconductivity with magnetically mediated 
pairing mechanism. Therefore, the fact that superconductivity appeared in the neighbor of 
the CO region is not only rare for organic conductors, but also important for possibilities of 
the superconductivity mediated by the charge fluctuation. Recently, SC state appeared in the 
neighboring region of the CO state is also found in /?- vanadium bronze, but no theoretical 
studies have been carried out. 

From the above background, we study /3-(DMcET)2PF6 in the following two point of 
views: i) how the checkerboard-type charge ordering is stabilized in the weakly dimerized 
organic system, ii) how the SC region appears in the neighbor of the CO region from the 
metallic state. 

The rest of this paper is organized as follows. In Sec. 2, we formulate susceptibilities 
and Eliashberg equations based on the RPA for the weakly dimerized 2D system. In Sec. 3, 
we determine appropriate parameter set of the model, and present results of the numerical 
analysis. Finally, summary and discussion are given in Sec. 4. 

2. Formulation 



In this section, we perform the analysis for the Hubbard-type model based on the random 
phase approximation following Kobayashi et alJ'^ Throughout this paper we set /cb = ^ = 1- 
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(b) Checkerboard (c) Stripe 




Fig. 1. (a) Model of the conduction layer for /3-(DMeET)2PF6 salt. The unit cell contains two 
molecules labeled by 1,2 with transfer integrals: tri, tr2, tqi, tq2 and tp. We also introduce cor- 
responding nearest neighbor Coulomb interactions: Ki, l^gi, 1^2 and Vp. Two candidates of 
the charge ordered states in the present analysis, (b) the checkerboard and (c) the stripe type CO 
states, respectively. The black and gray molecules are in charge-rich and poor states. 



2.1 Model 

We apply the extended Hubbard model on the two-dimensional (2D) lattice (see Fig. 1 
(a)), 

Tl = Tit + Tiu + 7iv , (1) 
E E(*-;^Wa'.S73<. + H-C-), (2) 

(ia;j/3) cr 

Ti-u = ''^Uaniainial, (3) 

ia 

Ti-V = E E ^ia;jl3''^iaa-njp^', (4) 
{ia\jl3) o-cr' 

where i,j (s I,-- - ,Nl) denote the lattice points of the 2D square lattice, and a,/? (g 1,2) 
specify molecules in the unit cell. c|q,^ {ciaa) is the creation (annihilation) operator for an 
electron with spin a (=t,i). {ia;jfi) represents a bond pair between the nearest neighbor 
sites, tia-j/3 denotes the transfer energy between sites {i,a) and (j, /5). TCfj and Tiy denote 
repulsive interactions where Ua and Va/s are the coupling for the on-site and those for the 
nearest neighbor sites, respectively. By using the Fourier transformation, 

Ciaa = -y= e 'Cj^, , (5) 

eq. (1) is rewritten as 



kk'qa 

2iVL ^ ^»l3^^)'^k+q,a,c7'^k,a,c7'^k' -q,l3,a''^k' ,P,a' (^) 
kk'qaficra' 
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Fig. 2. Energy band spectra and ^2 of eq. (13) and the Fermi surface. 



where the parameters in eq. (6) are taken as, en = £22, £12 = £21 ' ^11 = ^22, Vu — ^2*1- Then 
the transfer energies are 



£ii{k) = 2tp cos kx, 

£12 (fe) = trie-^'^y + tr2 + iqie''^- + tq2e-'^''^+''^\ 
The long range Coulomb interactions become 
Vn{q) = 2VpCOsqx, 

Vuiq) = Vrie-'iy + Vr2 + Vgie''!^ + Fq2e"'^«^+«^^ 
Applying the mean field approximation, the Hamiltonian is given by 

kajScr 



(7) 
(8) 

(9) 
(10) 

(11) 



(12) 



where a denotes the opposite spin of o. In this approximation, wc ignore the Fock term 
((cl^C g) ~ 0) and assume the uniform ground state: {uaa) = 3/4. The Hamiltonian is diago- 



nalized as 



(13) 



where Ci(fc) > E,2{k) ( sec Fig. 2 ). da-yaik) is the element of unitary matrix obtained by 
diagonalization of eq. (11). Then, {riaa) is given as 



^"'^"^ " ^ exp[(e,.(fe)-/.)/r] + r 



(14) 



where /x is the chemical potential determined by the condition of 3/4-filling: \ 5^q,ct('^ckt) = §• 
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2.2 Charge and spin susceptibilities 

We define matrices of susceptibility as 

fl/T 



-/ 



dr e'^"" {riaa (q) (r)n^^. (-q) (0)) , (15) 

where ui is the Matsubara frequency for bosonsJ r is the imaginary time. The density operator 
riaAl) is given by 

riaAq) = Yl 4,a,aCfe+q,a,a- (16) 
k 

The charge susceptibihty X'^ and the spin susceptibiUty for the easy axis are given by 

= 2 (^TT + ^iT + ^Ti + ^ll)a/3, (17) 

= 2 (^TT ~ ^iT - ^Ti + (18) 
Applying the random phase approximation and ignoring the Fock term, we obtain 

X,, = + -x(o)y(2)(/ + x(o)yW)-i 

xX(0)y(2)]-il(0)^ (19) 

X-^ = _(/ + x(o)y(i))-ixWy(2) [/ + xWy(i) 

-X(0)y(2)(7 + ^(0)y(i))-i^(0)^(2)]-i^(0)_ (20) 

and tX^) 

are given by 

vW = V + u^ -1/(2) =1/ (21) 

where V and C/ are defined as 

In this case, we treat mainly the nonmagnetic state, i.e., X^^ = X^, X-j-j = X^^. Then X'^ 
and are given by 

= {i + {U + 2F))-ix(o) , (23) 
X^ = X-f-f - 

= (/_X(0)[/)-ix(0), (24) 
where X*^ and X^ are hermitian matrices. The irreducible susceptibility X^^^ is 

= m ^S(^ + ^' i'^"* + ie„)G^°2(^' i^n) 



k.n 
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where e„ denotes the M 
Green function given by 



where e„ denotes the Matsubara frequency for fermions. G^^p{k,ioJn) is the single particle 



where ^jik) and da,j{k) are defined in eq. (13). 

Now, we introduce the linear combinations of the density operators as 

n±{q) = ^{ni{q)±n2{q)). (26) 

Then replacing the original density operators in eq. (15), and taking account the symmetry 
in the unit cell, the corresponding charge {v = c) and spin {v = s) susceptibilities are given 

by 

= X{'i±ReXr2. (27) 

The stability of the CO states in the present analysis (see Fig.l (b), (c)) is given by divergence 
of XI. 



2.3 Pairing interactions and Eliashberg equations 

The superconducting transition point is determined by the Eliashberg equation. The lin- 
earized Eliashberg equation for the singlet SC state is given by 

^ k',n',a',l3' 

xG^^l{k',ien')G%{-k', -i6„0S^V;/J'^(fe'), (28) 
with the paring interaction, 

= [7 + t> + ^UX^U - ^(f7 + 2V)X''{U + 2V), (29) 

where and X^ are the susceptibility matrices given in eqs.(23) and (24). In eq. (28), 
S^^.g-(A;) is the anomalous self-energy with the space inversion symmetry S^^.^-(fc) = 
S^g..Q^(— fc). Here we neglect the dependence of the Matsubara frequency in S^^.^-(fc), since 
such a treatment is valid for the weak coupling case.^"*^ In eq. (28), As = 1 corresponds to the 
superconducting transition point. 

Similarly, the linearized Eliashberg equation for the triplet SC state is obtained as 

ATS^,.^,(fc) = E Pla-M^-k') 

^ k',n',a',l3' 

X Gi°i (fe', ien')Gf^, i-k', -ien')K'a-Mk') , (30) 



6/15 



J. Phys. Soc. Jpn. 



Full Paper 



with the paring interaction, 




(31) 



In order to confirm the relevant part of the effective interaction for the singlet state, we 
divide the paring interaction into the following two parts. 



In the present linearized Eliashberg equations, the amplitude of the anomalous energies 
can be chosen arbitrary, so that we calculate the quasi-particle bands just below the SC 
transition temperature. The quasi-particle bands is obtained by the following Hamiltonian, 



where ai = a {ai = a) for the singlet (triplet) SC state. 

3. Numerical Results 

Since there arc a lot of parameters in the present model, we select a set of parameters in the 
following way. First, we fix the transfer integrals calculated by the extended Hiickcl method 
based upon the crystal structure analysis: tri = —0.0824, t,.2 = — 0.226, tp = 0.0475, tgi = 
—0.0438 and tq2 = —0.115 [eV].^^ In this case, we do not consider explicit pressure dependence 
of the parameters. Next, we determine the value of the on-site Coulomb energy U which 
becomes the ground state nonmagnetic. Then we determine the value of the nearest neighbor 
interactions V which stabilize the checkerboard type CO state. We use the unit for the energies 
as [eV] . Finally, we determine the phase diagram of this model. 

3.1 Nonmagnetic state 

To obtain the phase diagram in the (U, T) plane (Fig. 3(a)), we calculate which is larger 
than X^_. Figure 3(b) shows the momentum dependence of the spin susceptibilities X^, at 
T = 0.1. When the system is dimerized, the exchange coupling J = ^ becomes Jx > Jxy > Jy, 
where \tx\ = \tp - tqi/2\ = 0.069 > \txy\ = \tq2\/2 = 0.056 > \ty\ = |t,.i/2| = 0.0412.2 ^ 
result, X^ becomes the largest at (7r,0). /3-(DMeET)2PF6 is nonmagnetic. Therefore, we fix 
the on-site Coulomb interaction U at 0.4 [eV] in the following calculation where the spin 
density wave (SDW) does not appear. 

3.2 Charge ordered states 

We investigate the combination of the nearest neighbor Coulomb interactions V, where 
the checkerboard type charge order is stabilized. First, we examine the charge order patterns 




(32) 



(33) 



= HuF - XI ^iaa^kaa 



(34) 
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0.5 0.6 0.7 0.8 0.9 1 (0,0) (n.O) (jr^) (0,0) (-Jt.iT) (o.jr) (0.0) 

U 9 

Fig. 3. (a) Phase diagram in the {U,T) plane. The soUd hne correspond to the phase boundary 
between the metalhc and the SDW states, (b) Spin susceptibihties X^{q) for U = 0.4 and 0.5 at 
T = 0.01. 

(a) 
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( 

Fig. 4. (a) Phase diagram on (U, ^2) plane at T = 0.1. The dashed and the solid lines correspond 
to the CO and the SDW instabilities, respectively, (b) Charge susceptibilities Xliq) for {Vr2,l) = 
(0.35, 0), (0.294, 3) at T = 0.1, U = 0.4. 

by using the point charge approximation: Vi (x ^. Figure 4(a) shows the phase diagram on 
the (U, Vr2) plane at T = 0.1. becomes larger in the intermediate region between (0,0) 
and (— 7r,7r) as shown in Fig. 4(b), therefore we find that this result does not accord with the 
experimental fact. Taking account of the effective molecular size /, the coulomb interaction 
in the dimer is given as Vr2 oc (r^2 + ■ It is reported that / ~ 3 A is reasonable for 

the ET salt.^^ Therefore we fix / at 3 A and apply the point charge approximation for the 
other Vs. Figure 4(a) shows (C/, V^2)-pbase diagram at T = 0.1. In the case of / = 3 A, the 
CO instability occurs in the smaller Vr2 than that of / = A. This indicates that the charge 
fluctuation becomes larger by taking account of the molecular size, because the difference 
between Vr2 and another V is smaller. However, the critical point of Xt{q) does not change 
and we cannot obtain the checkerboard type charge ordering. 

In order to stabilize the checkerboard type CO state, it is suitable to satisfy the follow- 
ing condition: (K-25^7K?2) > (K-i)^gi)- For simplicity, we assume only two sets of nearest 
neighbor interactions Vi{= Vp = Vq2) and V2(= Vri = Vqi). As shown in Fig. 1(b), Vi is 
the interaction to induce the charge disproportion with the charge density in the unit cell: 
nn 7^ ni2- V2 is the interaction to favor the uniform charge density in the unit cell: rin = ni2. 
We set Vr2 = si X U, Vi = Vp = Vg2 = S2 x Vr2, V2 = Vri = Vgi = S3 X Vr2- For nonmagnetic 
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Fig. 5. (a) Phase diagram on {32,83) plane at T = 0.1, U ^ 0.4,^2 = 0.8U. The solid 
line corresponds to the CO instability, (b) Charge susceptibilities X'i{q) for (§2,33) = 
A(0.3, 0.8), B(0.4, 0.84), C(0.6, 0.78), D(0.7, 0.41) and E(0.6, 0.25) at T = 0.1, U = 0.4, Vr2 = O.m. 
The arrow shows the shift of the peak of X'l (q) from (0, 0) to (tt, tt) 





states, we fix [/ = 0.4 eV and Vr2 = 0.8 x U as shown in Fig. 3(a). 

Figure 5 shows the phase diagram on the {s2, S3) plane. As shown in Fig. 5(b), a peak of 
X'L appears at (0,0) in the stripe type CO phase (Fig. 1(c)), and (vr, vr) in the checkerboard 
type CO phase (Fig. 1(b)). The wave number for the peak of changes continuously in 
the intermediate region between the stripe and the checkerboard type CO regimes. In the 
parameter region of the point charge approximation, S2 is nearly equal to S3, but S2 should be 
sufficiently larger than S3 to realize the checkerboard type CO state. When we take account 
of the correspondence with the real material, it is suitable that the difference between S2 and 
S3 is small. Therefore, we set S3 = 0.3. 

Next, we examine the CO state by varying the temperature. The obtained phase diagram 
is shown in Fig. 6(a). There is a reentrant charge order transition as a function of temperature 
similar to the result of the extended Hubbard model in the 2D square lattice.'^'' We consider 
that this reentrant behavior is due to the competition between the effect of V{q) and X^^\q). 
At high temperatures, the peak of X^ is near (vr, vr). This originates from the momentum 
dependence of V{q). With decreasing temperatures, the momentum dependence of X^^\q) 
becomes large. Thus the peak of X^{q) changes from (vr,vr) as shown in Fig. 6 (b). 
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Fig. 6. (a) Phase diagram for the CO state on {s2,T) plane at U = 0.4, K2 = 0.8U, V2 = 0.3^2- The 
sohd Hne corresponds to the CO instabihty. (b) Charge susceptibilities X'i{q) for T — (A) 0.1, 
(B) 0.05 and (C) 0.01 at U = 0.4, Vr2 = 0.8f/, Vi = 0.64^-2 and V2 = 0.3^-2. 
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Fig. 7. Momentum dependence of the pairing interactions for the (a) charge [eq. (32)] and the (b) 
spin parts [eq. (33)], P^f{T), at U ^ 0.4, Vr2 = 0.8[/ and V2 = 0.3K-2- We set 53 = 0.64 at T = 
0.1 and S2 = 0.613 at T = 0.01, respectively. 



3.3 Pairing interactions and superconducting states 

First, we examine the effect of P^. Figure 7(a) shows the momentum dependence of P^. 
At {s2,T) = (0.64,0.1), Pfi shows a negative peak at (vr,7r), whereas P12 shows a positive 
peak. This result is consistent with the appearance of the checkerboard-type charge ordering. 
The effective repulsive interaction is caused by the Coulomb interactions where Vr2 and Vi 
are larger than V2. On the other hand, the effective attractive interaction arises by other 
Coulomb interactions V2. At {s2,T) = (0.64,0.01), the peak of P'^ is suppressed and changes 
from (7r,7r). We think it originates from the reentrant behavior of the CO. 

The momentum dependence of P^ is determined by that of X^^\ because U does not 
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Fig. 8. Phase diagram on (.S2, T) plane at U = 0.4, Vr2 = 0-8U, V2 = 0.3K-2- Tlic solid line corresponds 
to the CO instability. Note that the superconducting (SC) region is overestimated in the present 
analysis for the high-temperature region. The inset shows the detail close to T = 0.06 where the 
triplet SC state competes with the singlet SC state. 




Fig. 9. The momentum dependence of the element of anomalous self-energy matrix En of the singlet 
SC state for (a) T=0.01, S2=0.613 and (b) T=0.1, S2=0.64. 




Fig. 10. The momentum dependence of the element of anomalous self-energy matrix ReSi2 of the 
singlet SC state for (a) T=0.01, S2=0.613 and (b) T=0.1, S2=0.64. 

depend on the momentum. As shown in Fig. 7(b), the momentum dependence of increases 
with lowering temperatures. It is notable that is much larger than at T = 0.1 meaning 
that the SC state originates form the charge fluctuation. However, at T = 0.01, both and 
P^ are the same order. Therefore we expect that the SC state in this region is caused by both 
the charge and the spin fluctuations. 

In order to estimate the onset temperature for the SC states, we evaluate the linearized 
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Fig. 11. The momentum dependence of the element of anomalous self-energy matrix ImEi2 of the 
singlet SC state for (a) T=0.01, S2=0.613 and (b) T=0.1, S2=0.64. 




(b) 0.005 



■0.005 




Fig. 12. The momentum dependence of the clement of anomalous self-energy matrix En of the triplet 
SC state for (a) T=0.01, S2=0.636 and (b) r=0.1, S2=0.64. 




Fig. 13. The momentum dependence of the element of anomalous self-energy matrix (a) Re S12 and 
(b) Im E12 of the triplet SC state for T=0.1, S2=0.64 . 



Eliashberg equations eqs. (28), (30). Then we obtain the phase diagram on the (52,^") plane 
as shown in Fig. 8, the singlet SC state competes with the triplet SC state for 0.05 < T, but 
the singlet SC transition occurs first for T < 0.05. This means that an incommensurate CO 
phase is developed as lowering temperatures. 

As shown in Figs. 9 and 10, the anomalous self-energy of the singlet SC state strongly 
depends on the temperature. For the singlet SC state at T = 0.1, Re S12 is larger than En. 
This indicates that Re S12 contributes mainly to the SC state at T = 0.1. On the other 
hand. Re S12 is the same order of En. This indicates that both En and Re E12 contribute to 
the SC state at T = 0.01. As mentioned above, the spin fluctuation increases with lowering 
temperatures (Fig. 7) and is the same order of P'^ at T = 0.01. The pairing interactions 
do not have a special peak, since the anomalous self-energies are broad. 

The anomalous self-energy of the triplet SC does not depend much on the temperature as 
shown in Fig. 12. Re smaller than Pf^. It follows from eq. (31) that the total pairing 
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(b) 





?y - 





Fig. 14. The gap of the singlet SC state from the quasi-particle band of the noninteracting case at (a) 
T = 0.01 with S2 = 0.64 and (b) T = 0.1 with S2 = 0.613. The solid and dashed line correspond 
to the Fermi line and the gap, respectively. 




Fig. 15. The gap of the triplet SC state from the quasi-particle band of the noninteracting case at 
T = 0.1 with S2 = 0.613. The solid and dashed line correspond to the Fermi line and the gap, 
respectively. 



interaction of the triplet SC state is given by P'^ = P'^ — ^{P^ — U). This indicates that the 
triplet SC state is less affected by the spin fluctuation than the singlet SC state. Therefore, we 
consider that the triplet SC state is mainly induced by the charge fluctuation. It is reasonable 
that the triplet SC state appears near the CO state. 

Finally, we determine the symmetry of the order parameters by calculating the quasi- 
particle bands and count the numbers of zero points of gaps at the Fermi surface. We analyze 
eq. (34) as A ^ 1: we set A = 10~^ in this calculation, since the SC gap is regarded as 
infinitesimal. It is shown that the singlet SC state has 4 nodes as seen from Fig. 14. On the 
other hand, the triplet SC state has 2 nodes as shown in Fig. 15. In this result, however, we 
found a pathological behavior of the triplet gap which does not vanish at the center of the 
momentum space k = (0, 0), even though it satisfies the fermion antisymmetry. In fact, in the 
present analysis, we cannot classify the singlet and triplet gaps as the parity in the momentum 
space. This reason is explained as follows: The symmetry of the singlet (triplet) state is given 
as S^^(fe) = ibS^^(— fe), so that S^^(fc) = ibS^^(— fc) is not satisfied in general for a ^ (3. 
Therefore, the singlet (triplet) gap may include the parity-odd (-even) off-diagonal componets 
as shown in Fig. 11 (Fig. 13). In spite of this extra freedoms, the number of nodes coincides 
with the usual singlet- and triplet-gap properties. The explanation of this fact remains a future 
problem. 
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4. Summary and Discussion 

Wc have theoretically examined the CO and the SC states of /9-(DMeET)2PF6. We have 
analyzed the extended Hubbard model including the nearest Coulomb interactions in the 
weakly diniereized lattice. Applying the RPA, we choose the suitable parameters realizing the 
checkerboard type CO state under the nonmagnetic state. The CO phase shows reentrant 
behavior as a function of the temperature. We consider that this behavior originates from the 
competition between the wave vector dependence of V and the nesting effect of the charge 
susceptibility. 

Using the pairing interactions induced by the charge and spin fluctuations in terms of the 
RPA, we have estimated the onset temperature of the SC state near the CO instability by 
the linearized Eliashberg equations. The triplet SC state is stabilized around the region of the 
CO instability which originates from the wave vector dependence of V{q). On the other hand, 
the singlet SC state is stabilized near the region of the CO instability which originates from 
the nesting vector of the Fermi surface. 

The determination of parameters of the present model has been carried out to realize the 
checkerboard type CO state under the nonmagnetic state, comparing with the experimental 
result. However, values of these parameters are much different from those expected by the point 
charge approximation, even taking account of the effective molecular size. This suggests that 
the electron-phonon interaction due to the lattice distortion might be important to understand 
the checkerboard type CO state. 

In this paper, we have assumed that the CO state is an insulating state. There is a 
possibility, however, that the CO state is a metallic state due to existence of the hole and the 
electron pocket in the Fermi surface, and a SC state coexists with the CO state. This scenario 
was asserted to explain the superconductivity of a-(ET)2l3 salt under uniaxial pressure.'^'® 
Actually, the temperature and the pressure dependences of the resistivity in certain parameter 
regions seem to behave similarly in these two systems. "'^'^'^"'^ If such a metallic state is realized, 
the phase diagram of Fig. 8 may be modified so that a different type of SC state appears in 
the CO region, but the SC states obtained in the present analysis remain to be unchanged. 
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